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ABSTRACT 


With the strong CGRO/BATSE evidence that most 7 -ray bursts do not come from 
Galactic neutron stars, models involving the accretion of a comet onto a neutron star (NS) 
no longer appear to be strong contenders for explaining the majority of bursts. If this 
is the case, then it is worth asking whether the lack of an observed Galactic 7 -ray burst 
population provides a useful constraint on the number of comets and comet clouds in the 
Galaxy. Owing to the previously unrecognized structural weakness of cometary nuclei, 
we find the capture cross sections for comet-NS events to be much higher than previously 
published estimates, with tidal breakup at distances Uj % 4 x 10 10 cm from the NS \s 
a result, impacts of comets onto field NSs penetrating the Oort Clouds of other stars are 
found to dominate all other Galactic NS-comet capture rates by a factor of 100. This in 
turn predicts that if comet clouds are common, there should be a significant population of 
repeater sources with (i) a Galactic distribution, (ii) space-correlated repetition, and (ii) 
a wide range of peak luminosities and luminosity time histories. If all main sequence stars 
have Oort Clouds like our own, we predict ~4000 such repeater sources in the Milky Way 
at any time, each repeating on timescales of months to years. Based on estimates of the 
sensitivity of the CGRO/BATSE instrument and assuming isotropic 7 -ray beaming from 
such events, we estimate that a population of ~ 20 - 200 of these Galactic NS-Oort Cloud 
7 -ray repeater sources should be detectable by CGRO. In addition, if giant planet formation 
is common in the Galaxy, we estimate that the accretion of isolated comets injected to 
the interstellar medium by giant planet formation should produce an additional source of 
Galactic, non-repeating , events. Comparing these estimates to the 3-4 soft 7 -ray reeater 
sources detected by BATSE, one is forced to conclude that (1) comet impacts on NSs are 
inefficient at producing 7 -rays; or (2) the 7 -rays from such events are highly beamed; or 

(3) the fraction of stars in the Galaxy with Oort Clouds like our own is not higher than a 
few percent. 
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1. INTRODUCTION 


The first widely explored mechanism for producing astrophvsical 7 -rav bursts (here- 
after denoted GRBs) involved the accretion of comets onto a neutron star (Harwit & 
Salpeter 1973; Colgate & Petschek 1981; Van Buren 1981; Tremaine & Zytkow 1986; cf. 
the review by Higdon & Lingenfelter 1990). Most studies of this mechanism focused on 
events in which a parent neutron star (NS) accretes bound comets from a comet cloud 
surrounding its progenitor. However, two main variants on this theme have been reported. 
In the first (Tremaine & Zytkow 1986), the parent NS accretes comets from its own Oort 
Cloud (OC) when the cometary periastra are sufficiently small. This scenario is prob- 
lematic. however, since it requires that (i) high-mass stars form Oort Clouds, (ii) these 
stars retain their OCs after the NS-forming event, and (iii) a sufficient number of comets 
then pass within the NS capture radius to create the rate of observed GRBs. In a second 
scenario, Pineault (1990) offers a mechanism for removing concerns (ii) and (iii) by having 
a newly-formed NS capture comets from its own OC as it exits from the center of this 
cloud, after an asymmetric supernova (SN) explosion. 

A small fraction of the GRBs known as Soft Gamma Repeaters (SGRs) exhibit re- 
current activity (Atteia et al. 1987), softer spectra (Norris et al. 1991), and a possible 
association with supernova remnants and pulsars (Kouveliotou et al. 1993; Kulkarni Sz 
Frail 1993). Prompted in part by the evidence for a NS — SGR association, we explore 
here a third scenario, in which field NSs penetrate the OCs of other stars and accrete 
comets from these Extra-Solar Oort Clouds (ESOCs). This scenario was independently 
proposed by Pineault & Poisson (1989) and Mitrofanov (1990), and later expanded on by 
Mitrofanov & Sagdeev (1990) and Pineault (1992). 

Our contribution to this scenario is four-fold. ( 1 ) In §2 we re-evaluate the expected 
NS-ESOC comet accretion rates using more modern cometary and Oort Cloud parameters. 
Because the capture radius we derive is significantly larger than previously recognized, we 
find that the event rates produced by this mechanism dominate over captures of comets 
onto their parent NS by a factor ~ 10 ^. ( 2 ) Owing to the larger event rates, in §4 we 
find that there should be a significant number of Compton Gamma Rav Observatory 
(CGRO) detectable Galactic NS-comet accretion events each year, for which we describe 
the expected characteristics of NS accretion of both ESOC comets and interstellar comets. 
(3) Later in §4 we use the paucity of CGRO-detected SGRs to place a significant constraint 
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on the fraction of stars with cometary Oort Clouds like our own. (4) In §5, we examine our 
model to determine the number of GRBs that must be detected to constrain the population 
of interstellar comets. In estimating the capture cross sections and encounter statistics, 
we also consider the effects of thermal ablation and evaporation in the radiation field of 
the neutron star (Katz 1986). 


2. ENCOUNTER RATE ESTIMATES 

We begin by computing the encounter rate between field NSs of space density n na 
and ESOCs of density n,/ oc , where n* is the mean local space density of stars and f oc 
is the fraction of these stars with Oort Clouds. The local Galactic mass density (Bahcall 
1984) is (0.185 ± 0.02)M o pc -3 , of which approximately half can be accounted for by 
stars. If we adopt a mean stellar mass m = 0.5 A/©, then n, ss 0.2 pc -3 . For free 
NSs we adopt a local space density n n3 = 2 x 10~ 4 pc -3 , corresponding to a total NS 
population N na = 2.5 x 10 8 in the Galactic disk and implying a formation rate of 1 per 40 
yrs over 10 Gyr. We adopt a typical NS velocity u nj = (100 km s -1 )u 10 o. Our values for 
both n na and v na are consistent with recent orbital calculations (Blaes & Madau 1993) in 
which pulsars with a bimodal distribution of initial velocities (Narayan Sz Ostriker 1990) 
are injected into a Galactic disk/spheroid of total mass 8.78 x 10 10 A/©. For an (r,v) 
phase-mixed NS vertical distribution with a half- number scale height (- 1 / 2 ) = 280 pc and 
a root-mean-square height z rma = 710 pc, a mean stellar mass of 0.5 A/© implies that 
the total number of Galactic disk stars is JV* = 1.76 x 10 n , and the effective Galactic 
volume is \ gal = 9.5 x 10 11 pc 3 . The same assumptions give a velocity distribution with 
(u) = 81 km s -1 and v rma = 98 km s -1 , which is somewhat slower than the injected pulsar 
distribution. In § 6 we discuss possible changes in the encounter rates and accretion arising 
from a recent, upward revision (Lyne & Lorimer 1994) of the pulsar velocity distribution. 

To investigate our scenario quantitatively, we use a simple kinetic model to compute 
the NS - ESOC encounter rate over the Galactic disk. This yields 

Nns = n na (7rR^ c )u nj iV. = (10 yr -1 )/o C uioo-R6ooo • (1) 

Here we have adopted a typical ESOC core radius R oc = (6000 AU)R 600 o based on R oc % 
7500 AU for the dense inner core radius of our own OC (Weissman 1990), scaled to a mean 
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stellax mass of 0.5 M 0 . Although we could have added the extra parameter of specifying 
the radial distribution of comets within the ESOC (Pineault & Poisson 1990), we doubt 
whether that more sophisticated approach would result in any improvement in the viability 
of the scaling relations. Importantly, each NS - ESOC penetration will last a sufficiently 
brief time. 

Tns ~ -z- 1 — = (400 yrJ/Jeooowfoo > (2) 

that individual SGR sources should turn off after a few centuries. The intruder NS-ESOC 
mechanism therefore predicts that one of the 3-4 SGRs now known should decrease or 
cease its burst rate within 50-100 yr. 

Combining equations (1) and (2), we find that the total number of intruder NSs in 
ESOCs spread throughout the disk at any given time should be 

= N na T n s = Tl ns \ oc N + f oc ~ 4000/ oc i?g 00 g , (3) 


independent of the NS velocity, where V oc is a typical Oort Cloud volume. 

By comparison, consider the number of active sources in which a NS accretes comets 
during post-formation ejection through its own OC (Pineault 1990). For a typical exit 
velocity of 100 km s -1 , a NS should spend approximately (300 yr)R 600 o departing its own 
OC core. Because this is ~ 10 -7 of the 3-5 Gyr age of a typical NS, there should exist 
an average of 10“' / n3iOC iV na , or ~ 25 fns.oc active sources of this type in the disk, where 
fns.oc is the fraction of NS progenitors that have Oort Clouds. This ratio can be expressed 
more simply as 


C = 10 7 




(4) 


This simple comparison demonstrates that there should be ~ 10 2 times more intruder NSs 
inside the Oort Clouds of other stars than those in their own Oort Clouds at any given 
time. 

The rate of NS-comet encounters in OCs around parent NSs can only dominate the 
rate due to NS encounters in the ESOCs of other stars if either: (i) the fraction of NSs 
with OCs exceeds the average fraction of all stars with OCs by a factor 10 2 , or (ii) NS 
Oort Clouds have ~ 10 times larger volumes than field-star Oort Clouds. In our view, 
neither of these alternatives is likely (cf., also Tremaine 1993). Therefore, the intruder 
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NS rate of encounters almost certainly the dominant source of NS-comet collisions in the 
Galaxy. This is our first significant result. 

Next, we desire to predict the burst frequency of NS - ESOC sources. As a precursor 
to that, however, we must compute the NS - comet capture radius, Rq. We will do so 
using a new capture mechanism which relies on the tidal breakup of fragile icy comets (see 
§ 3). Based on tidal breakup or on the drag introduced by the motion of a conducting 
body through the NS magnetic field, the capture radii in previous models of NS - comet 
encounters were of order 3 x 10 7 cm (Van Buren 1981; Tremaine & Zytkow 1986; Pineault 

1990) . To investigate the tidal breakup of comets, Colgate & Petschek (1981) assumed a 
comet tensile strength s = 10 10 dyn cm -2 , characteristic of a Ni/Fe composition, rather 
than the more fragile “icy mudballs” in the standard comet model (McDonnell et al. 
1993). Later, Tremaine & Zytkow (1986) adopted tensile strengths between 10 6 (Wetherill 
& ReVelle 1982) and 10' dyn cm -2 (Hobbs 1974), which is appropriate for solid ice. 
However, evidence based on observations of Sun-grazing comets (Sekanina 1982; Hughes 

1991) points to strengths in the range of 10 3 < s < 10 5 dyn cm -2 ; more recently, the 
breakup of comet Shoemaker- Levy 9 around Jupiter (Sekanina 1993; Weaver et al. 1994) 
suggests that s ~ 10 3 dyn cm -2 (Scotti & Melosh 1993). In what follows, we conservatively 
adopt s = (10 5 dyn cm _2 ).S 5 . 

Comets in our Solar System have typical masses between 10 16 g and 10 18 g, and consist 
of comparable amounts of water ice and heavy elements (Weissman 1990). We consider 
a canonical comet of mass m c — (10 17 g)mi 7 , density p c = (1.0 g cm -3 )p 0 , and radius 
R c = (2.88 kmjm^ p 0 approaching a NS of mass M ns = 1.4 M®. The tidal breakup 
radius for such comets in the NS gravitational field is 

Rb = ( Po RiGM na /s) 1/3 zz (5.4 x 10 10 cm) m 2 7 9 p 1 ^ 9 . (5) 

If we consider the gravitational focussing of the comet’s orbit by the NS, a typical approach 
velocity v n3 , Rb corresponds to an impact capture parameter of 

Ro = (2GM ns R b ) 1 ^/v na « (4.5 x 10 11 cm) m| 7 /9 Voo • (6) 

Most evidence suggests that comets are a conglomerate of fragments with varying tensile 
strengths, including gravel, dust, and ice particles. The weakest fragments will be broken 
at the largest distances, followed by a cascade down to the size of particles of sufficiently 
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small size and high tensile strength (i.e., 10 7-9 dyn cm -2 ) to resist the tidal forces. We 
expect this cascade to occur over time scales comparable to the sound-wave crossing time 
of the solid fragments (~ seconds). Because the strength of the individual fragments is 
expected to be much weaker than ice, the tidal cascade should not cease until the fragments 
reach the size of strong (i.e.. molecular) grains of ice and dust. We show in § 3 that these 
grains will quickly evaporate and become ionized in the NS radiation field to create a 
plasma. 

Once the cascade to smaller fragments begins, the debris area rises substantially, so 
that significant radiation can be absorbed. The difficulty with this cascade from km-sized 
comets to cm-sized fragments is that R b oc (^/s) 1 / 3 . Thus, smaller fragments eventually 
become stable against tidal forces, but the factor of 100 range in tensile strength, s = 10 3— 5 
dyn cm only translates into a factor of 5 in breakup size. However, because the process 
of fragmentation is incompletely understood and may involve fragment collisions or forces 
of evaporative jetting, some unknown fraction of the fragments will probably continue on 
hyperbolic trajectories. 

Again adopting a kinetic model, we predict that the comet - NS collision rate in each 
of the 4000/ oc active Galactic NS - ESOC sources is 


N c = n c (TrR-)v n3 = (0.65 yr _1 )JV : 


13 


-1 


E? — 3 - , 

- ft 6000 U 100 m l7 


2 / 9 c ~ 1 / 3 . 1 / 9 


P 0 


( 7 ) 


where n c is the average number density of comets in the comet cloud. R 0 is the NS- 
comet capture radius, and A r i 3 is the number of comets per OC in units of 10 13 , the 
estimated population for our own OC (Weissman 1990; Mumma, Weissman. &: Stern 
1993). Multiplying equations (3) and (7), we see that the total event rate in the Galaxy is 
of order 2600 yr -1 , some fraction of which are near enough to be detectable (see § 4). 

Clearly, with a typical encounter timescale of a comet everv ~ 1 — 2 vrs, field NS 
penetrations of ESOCs will produce repetitive, spatially correlated accretion events. Be- 
cause the details of cometary breakup and accretion onto a NS are strongly dependent 
on the structural properties of the nucleus, the NS mass, and the encounter geometry 
whether or not some debris is trapped in a temporary orbit in the NS magneto- 
sphere), we expect a wide range of accretion luminosity time histories to result, ranging 
from milliseconds to minutes. 
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3. COMET - NEUTRON STAR COLLISION CROSS SECTIONS 


As noted in §2, previous NS-comet tidal disruption studies adopted tensile strengths 
that are 3 - 5 orders of magnitude larger than those in the standard '‘icy mudball” comets 
(Mumma et al. 1993; Sekanina 1993) seen around our Sun. The result of a NS-comet 
encounter is a debris cloud of mass m c , initially moving along the original hvperbolic 
approach trajectory at breakup, of velocity 


Vb 


RftVns 

Rb 


(830 km s" 1 ) m“ l/9 s\ /6 p~ 1/is 


( 8 ) 


For impact parameters less than R 0 , the comet will penetrate inside R 6 , where it will 
break into a number of smaller pieces, each with a random velocity about v b provided by 
differential tidal acceleration and heating-related jetting forces. Charged particles will be 
affected by the NS magnetosphere. In this section we calculate the tidal capture radius for 
icy mudball comets. We then discuss the fate of this tidal debris in the strong gravitational 
and magnetic fields, and the intense X-ray/UV radiation environments of NSs. 


3.1 Ablation and Evaporation 


We now show that during the inbound portion of the encounter, of duration 

U nc = (nRb/vb) = (2000 s) m\ /3 s~ l/2 p \ /6 , (9) 


significant evaporation of the tidally-produced fragments occurs in the NS radiation field. 
Once the breakup process begins, the surface area exposed to NS heating increases. The 
flux of NS radiation is set by its surface area and surface temperature: the latter is de- 
termined by accretion from the interstellar medium to be least 3 x 10 5 Iv and perhaps as 
high as 10 6 K (Wang & Sutherland 1994). If we parameterize T ns = (10 5 ' 5 KjTs.s and 
L ns = (7 x 10 30 ergs s -1 )T 5 4 5 , the blackbody radiative flux at a distance R = R b from a 
10 km radius NS is 


F — 

■L ns — 


^ = (2 X 10 8 ergs cm" 2 s' 1 ) ^ . 


( 10 ) 


For large fragments, in which thermal conduction across their diameter can be neglected, 
the equilibrium temperature is 


T c = 


R n 


R b 


1/2 


1 - A 

4e 


1/4 


T aa (980 K) T 5 . 5 m- 1/9 4 /6 p- 1/18 . 


( 11 ) 
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for a typical comet albedo A « 0.04 and emissivity e « 0.9 at R = R h . At this high 
temperature, significant sublimation of cometarv ice will occur, thereby limiting the surface 
temperature. By equating the latent heat flux of sublimation to the radiative heating flux, 
for values given in equation (10) with scaling factors set to unity, we estimate that the 
actual temperature is ~ 240 K, with a mass loss rate per unit area ~ 0.002 g cm" 2 s -1 . 
This is sufficient, over the encounter time, to erode several centimeters from the surfaces 
of all comet fragments. 

To estimate the fragment size that can sublimate completely, we first calculate the 
^fluence” (energy A E per projected area A) delivered to a single comet over the total 
encounter. For an orbit with specific angular momentum Rov ns , this is given by 


A E 
A 


Lns 

2i? 0 ^r 


= (8 x 10 n ergs cm 2 ) T * 5 


m 


- 1 / 9 1 / 6 - 1/18 


17 


P 0 


(12) 


which is sufficient (Katz 1986) to ablate cometary surfaces to a depth 


A R 


(A E)/A 
4 p c h{T) 


V np 4 - 1/9 1/6 - 1/18 

= (7 cm) T 5 _ 5 m 17 s 5 ' p 0 


(13) 


where h(T) « 3 x 10 10 ergs g _1 is the steam enthalpy at T « 240 K. Radiative cooling can 
be shown to be negligible under these conditions. 

The ablation depth is insensitive to the comet parameters, but is extremely sensitive to 
the NS radiation temperature, which is set by the rate of interstellar accretion, els discussed 
further in § 3.2. The NS surface temperature of 10 5,5 K corresponds to a mass accretion 
rate M = 8 x 10 10 g s -1 onto the entire NS surface (10 km radius) at 10% efficiency. Such 
accretion rates are entirely consistent with the detailed modeling of NS accretion from the 
ISM (Blaes & Madau 1993), although large differences in the accretion rate and X-ray 
luminosity can arise for different phases of the ISM and from polar-cap accretion. The 
population of SGRs might, for example, arise from NSs accreting from regions with higher 
than average density. Old NSs might also be expected to have somewhat lower magnetic 
fields, suggesting that the higher X-ray luminosities and radiation temperatures associated 
with polar-cap accretion are not present. 

In the NS radiation field assumed here, fragments of radii ^ 0.1 m will be completely 
vaporized. Fragments as large as several meters could evaporate only around NSs at 
10 6 K. The fate of the larger chunks of the comet cascade is uncertain and probably 
depends on processes that occur during the tidal breakup, such as ablation, jetting, and 
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fragment collisions. The limiting fragment size may drop from 1 km to 1 m in proportion 
to (i? 2 /^) 1 / 3 . These processes may be able to bridge the gap between the initial km-sized 
fragments of the comet and the sub-meter sized fragments than can evaporate in the NS 
radiation field. Thus, a significant fraction of the comet mass may be able to attach to the 
NS magnetic field lines and accrete to the NS. From equations (12) and (13), we note that 
significant fragment ablation (to depths A R on the order of meters) requires fluences r 10 13 
ergs cm -2 corresponding to a NS surface temperature 10 5 ' 8 K at an impact parameter Ro 
given by equation (6). Alternatively, we can re-write the impact parameter necessary for 
a given Ail as 

Ro = (3 x 10 cm) 5 v 100 p 0 ^ ^ ' (14) 

Thus, including gravitational focussing, we have demonstrated that comets of tensile 
strength s = 10 5 dyn cm -2 are tidally disrupted for impact parameters i? 0 < 4.5 x 10 11 
cm. For impact parameters R 0 ranging from 4 x 10 11 cm inward to 3 x 10 10 cm, increasing 
amounts of the tidal fragments will be ablated. In the intense radiation field near the 
NS, the molecules produced by the sublimation of the cometary ice and dust fragments 
will be quickly dissociated and ionized. The photoionization time of H 2 O and its primary 
daughter products OH and H at any distance R from the NS is quite short, 

f ' S5(3s) (^) 2 (15) 

where R6,o = 4.5 x 10 10 cm. as given in Equation (6). The resulting plasma should consist 
of ions that can attach to the NS magnetic field lines and accrete to the NS surface, since 
the effects of NS radiation pressure can be shown to be small. 

3.2 Magnetospheric Capture 

The processes governing the magnetospheric penetration of this debris cloud and its 
subsequent accretion onto the NS now deserves discussion. To gain some insight into which 
processes dominate, we estimate criteria for: (1) the debris cloud avoiding magnetospheric 
deflection; and (2) the plasma attaching to the magnetic field lines. We assume that 
the plasma from the comet is quickly heated to temperatures ~ 10 4 K by the UV/X-ray 
radiation from the NS. For a silicon/oxygen plasma, this temperature corresponds to a 
mean thermal expansion speed of u exp « (3 km s -1 )t> 3 ; for H the thermal expansion speed 
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is ~ 10 km s . Over the characteristic duration of the encounter. t e nc , this plasma cloud 
will expand to fill a sphere of radius R exp % t 7 exp t enc . 

To set the first criterion for whether this diluted debris cloud is deflected from the 
magnetosphere or instead penetrates to the NS surface, we compare the ratio of the ram 


pressure of the plasma debris cloud. pV 2 . to the magnetic energy density, (B 2 /Sir), at 
the breakup radius R b . For the NS magnetosphere, we assume a dipole magnetic field, 
B(R) = B 0 (R ns / R) 3 -, typical of an old NS. Current studies show that these fields may be 
high. Whereas the majority of early pulsar studies assumed that the initial large ( IO 125 G) 
magnetic field decays exponentially on a time scale of order 10 7± ° 2 yr (Gunn & Ostriker 
1970: Lyne, Manchester, Taylor 1985), more recent investigations (Narayan & Ostriker 
1990) show that the form of the field decay (exponential or power-law) is uncertain. There 
also have been suggestions for both a long-lived residual component B < 10 10 G (Kulkarni 
1986) and for little or no field decay in isolated NSs over 10 8 yr (Bhattacharya et al. 1992; 
Chanmugam 1992). We adopt a surface field parametrization Bq = (10 10 G)Z?io, realizing 
that the factor B io may range between 10 -2 and 10 2 . 

For the ram pressure of the comet debris cloud, we assume that the plasma de- 
bris cloud moves at bulk velocity v = v b and that the debris cloud density is p b = 
m c/[(l ,r /3)(^exp^enc)' ! ]- The much lower random velocity of the solid fragments is pro- 
duced by fluctuations in the tidal acceleration. Since most of the fragments are vaporized, 
we can use the fluid approximation to evaluate the ratio of ram pressure to magnetic 
pressure at R b , 


B 2 /8tt 


= (5 x 


i a 15 \ iu/y — i/u i/ioj-j 

10 ) 771 | 7 pQ V‘ 


1/18 0—2 — 3 

10 U 3 


(16) 


Because the above ratio is much greater than unity for any likely NS magnetic field, comet 
mass, density, and tensile strength, the plasma debris should easily penetrate the NS mag- 
netosphere. Of course, some of the larger comet fragments may escape full vaporization, 
and therefore avoid capture by the NS magnetosphere. 

We now investigate the second criterion, which governs whether the accreted plasma 
attaches to the magnetic field lines of the NS. The Larmor radius of a debris fragment of 
mass Nm v , charge Ze, and gyrovelocity v b arising from the original comet of mass m c is 
given by 


Rl 


[Nrn^cvhR^ 

( ZeBoR l s ) 


«(1 x 10 8 cm) 



m 


5/9 —5/6 
17 ^5 


5/18 

P 0 


R iq 1 


(17) 


This gyro-radius is considerably less than R b for fields B 0 > IO 10 G. and can be reduced 
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for values of N/Z appropriate to a vaporized fragment. Therefore, the accreting plasma 
will follow the magnetic field lines down to the NS surface, once it penetrates the magne- 
tosphere. 

We evaluate the radius. Rm , of the magnetosphere first for a NS accreting at rest 
from a static ambient ISM with isothermal sound speed c a = (10 km s _1 )c 10 , and then for 
a NS moving at velocity v na c 3 . For an isolated, low-velocity NS, we assume isothermal, 
spherical accretion flow and equate the ram pressure, pv 2 = {p a c 2 ){r j R a )~ 5 G ^ of the 
accreting matter inside the sonic radius, R a = (GM na /2c 2 a ), to the magnetic pressure, 
{B'l /S^){r j R na )~ b . The density at the sonic radius, p s = e 3 / 2 p 00 , can be related to 
the mass infall rate, M = (e 3 / 2 /4)(47r)p 00 (GM na ) 2 /c^, and to the NS luminosity, L na = 
M(r} a c 2 ), with 77 a zz 0.1. We then find a magnetospheric radius independent of c , , 


Rm = 


L na (2GMns) l/2 \ 


2/7 


= (1.2 X 10 9 cm)flV’:r.-«/T , 


(18) 


where we have expressed L na in terms of the NS effective temperature. For the more likely 
case where v na c a , one substitutes v na for c a (Hoyle Lyttleton 1941). The formulae 
for accretion radius and Rm may be written in terms of the NS luminosity, similar to eq. 
(18), or in terms of ambient density and NS velocity, 


R* 


GM„ 


= (1.9 x 10 


12 


cm ) u 100 


(19) 


and 


Rm = 


b WA 


2/7 


= (4 x 10 £ 


cm)B 4 ^ 7 ’' 6 ^ 


'10 


u ioo 


- 2/7 


(20) 


Evidently the tidal breakup and fragment vaporization occur at distances ~ 10 10-n cm 
from the NS, which is larger than the magnetospheric radii for old NSs accreting at a 
standard rate from the ISM. Because the Laxmor radii of the evaporated fragments are 
much less than both Rb and Rm , it is likely that the ablated and ionized debris attaches to 
the field lines of the NS. Although the comet debris cascade will be complex, it is possible 
that significant amounts of the cometary dust and ice axe broken up, sublimated and 
ionized so that they penetrate the NS magnetosphere and accrete onto the NS to produce 
a 7-ray burst as previously explored (e.g., Colgate & Petschek 1981). If the capture 
radius is i?o ~ 4 x 10 11 cm, as argued in equation (6), then the encounter frequency is 
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interestingly large (~ 1 yr 1 ). If the effective Ro is a factor of 4-5 smaller, owing to 
incomplete evaporation of the comet fragments, the encounter rate drops to a value 

= (0.033 yr -1 )iVi 3^000 ( 1Q n° cm ) ( 21 ) 

In that case, the repetition frequency of bursts might be dominated by bound comet debris 
returning to the NS on later orbital passages. 


4. COMET - NS IMPACT DETECTABILITY AND CHARACTERISTICS 

We now consider the detectability of cometary accretion events onto NSs. For a comet 
of mass m c = (10 17 g)mir, the energy released is a significant fraction of the comet’s rest 
mass, and may be written 

E = rjm c c 2 (9 x 10 36 ergs) m 17 (^-) • (22) 

Here, 77 is the net 7-rav conversion efficiency of the collision, normalized to the comet 
rest-mass energy. Since the gravitational potential energy at the NS surface is about 20% 
of the rest-mass energy, 77 = 0.1 corresponds to conversion of half the kinetic energy of 
infall into 7-rays. Values a factor of 5 lower have also been used (Tremaine Sz Zytkow 
1986). Since the CGRO/BATSE 7-ray burst detector has sensitivity to burst fluxes down 
to Fmm = 3x 10“ 8 ergs cm" 2 s" 1 (Schaefer 1993a; Meegan et al. 1992), a comet accretion 
event onto a NS could be detected to a distance of 

j j 2 

' Uwr) «a6kpc)m^(^ I ) 1/ hiT)-/ 2 , (23) 

where 6T is the burst duration timescale in seconds. For 6T = 0.1 s (Kouveliotou et al. 
1993; Higdon & Lingenfelter 1990), we obtain D max = 5 kpc. Although our poor knowledge 
of 77 m ak es these estimates highly uncertain, both the 1.6 kpc and 5 kpc detection distance 
scales are small compared to the radius of the Galactic disk, but large enough to reach 
neighboring spiral arms. Since D m ax is comparable to the vertical scale height of the “thick 
stellar disk” and the 710 pc thickness of the phase-mixed NS population (Blaes &: Madau 
1993), NS-comet gamma ray bursts produced by field NS intruding into ESOCs should 
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exhibit a Galactic distribution. Since this is contradiction to the BATSE observations for 
the general GRB population, we can obtain a constraint on f oc for the SGR population. 

Assuming that the spatial distribution of NSs and OCs are each homogeneous over 
the Galactic disk, then our estimate of ~ 4000/o C active NS — OC sources in a disk of 
radius 12 kpc, detectable to 1.6 kpc, implies ~ 70 f oc BATSE-detectable, repeating NS - 
OC sources at any given time. For the more relevant = 0.1 s SGR burst timescales 
detected by BATSE, corresponding to a detection limit of 5 kpc, we predict ~ 700/ oc 
BATSE-detectable, repeating NS - OC sources at any given time. This is our second 
result. 

Only three (Hurley 1993) or four (Wang &: Lingenfelter 1993) SGRs have been identi- 
fied to date. Two of these sources can be associated with NSs: SGR 1806-20 is associated 
with a Galactic pulsar and supernova remnant (Kulkarni & Frail 1993), while SGR 0525-66 
has been associated with the LMC supernova remnant N 49 (Mazets et al. 1982; Helfand 
& Long 1979) and a compact X-ray source, possibly a pulsar (Rothschild et al. 1994). 
Recurrent activity over more than 20 yrs in the repeater SGR 1900+14 (Kouvelioutou 
et al. 1993) and repetitive soft 7 -ray bursts from GB 790107 (Livio & Taam 1987) lend 
further support to suggestions that SGRs are related to NSs. Thus, unless no 7 -rays are 
produced in comet accretion mechanism or unless the 7 -rays tire sufficiently beamed that 
few events axe detectable from Earth, the occurrence fraction of stars with Oort Clouds 
similar to our own is not high. As a result, we estimate that for D max = 1.6 kpc, no more 
than 3/70 « 4% of all stars have comet clouds similar in population and structure to our 
Oort Cloud. This is our third result. Because this limit depends crucially on 9 raax , this 
figure is only illustrative. 

If, in the future, the already-detected SGRs are found to reside beyond the ~ 1.6 
kpc upper limit on the (BATSE) NS-comet detection distance, then f oc can be even more 
tightly constrained to < 1 %. Since the formation of OCs like our own clearly depends on 
the presence and astrocentric distribution of giant planets in the distant solar system (cf., 
Weissman 1990; Tremaine 1993), our results suggest that such solar system architectures 
could be rare. In the same regard, it is worth noting that even if the impacting comets 
were as massive as 10 18 g, their collision with a NS would not likely be detected by BATSE 
in the LMC (e.g., the SGR tentatively identified with the N49 supernova remnant). On the 
other hand, detections of massive-comet captures by NSs in the nearest globular clusters, 
at 3 - 10 kpc, would be possible. 
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5. NS COLLISION RATES WITH INTERSTELLAR COMETS 


In addition to collisions between ESOCs and NSs. field NSs can also accrete free 
comets residing in the interstellar medium (ISM). The best estimates for our OC (Weissman 
1990: Tremaine 1993) predict that between 2 and 5 times the number of comets injected 
into bound orbits were ejected hyperbolically, with Jupiter and Saturn being the primary 
scattering centers for interstellar comets. It is therefore widely expected that the ISM 
contains a significant population of comets (Sekanina 1976; McGlynn & Chapman 1989; 
Stern 1990). 

If there are a > 1 interstellar comets (ISCs) for every OC comet in the Galaxy, 
then field NSs should encounter interstellar comets at a times the rate of OC comets. 
Adopting the same formalism as in Eqns. [3] and [7] above, we estimate that the ISC 
source should produce 2600 ocf oc NS-comet collisions yr -1 in the Galaxy. Thus, if comet 
clouds like our own are rare, there should exist a = 2 — 5 times as many isolated (i.e., 
non-repeating) Galactic GRB sources due to NS-interstellar comet collisions as there are 
NS - OC repeaters. However, for the same capture radius as in equation (6), we predict 
that for each NS, the mean time between interstellar comet collisions is in excess of 10 4 yr. 
Because dynamical mixing will cause the interstellar comets to be smoothly distributed in 
the ISM, the space density of interstellar comets will be ~ (2 x 10“ 5 )a times less than that 
of OC comets, and, unlike NS - ESOC collisions, collisions between NSs and ISCs will not 
be spatially clustered. This is an important distinction. 

We estimate that a Galactic population of between 90/ oc yr -1 and 230 f oc yr" 1 non- 
repeating sources produced by NS collisions with ISCs should be detectable in the BATSE 
dataset for a = 2 — 5 and a BATSE burst detection limit of 1.6 kpc. These sources should 
exceed the SGR bursts by a factor a = 2 — 5. As demonstrated by Smith &c Lamb (1993) 
and Syer & Saha (1994), one can hide up to 20% of the population of BATSE- detected 
GRBs distributed in the Galactic disk and remain consistent with current BATSE isotropy 
measurements. This would permit up to 80 possible Galactic NS - ISC events yr -1 to reside 
in the BATSE dataset, currently (June 1994) totaling over 1000. Therefore, as our fourth 
result, we conclude that the present BATSE dataset does not allow the population of ISCs 
to be meaningfully constrained. However, over the next several years, as the number of 
GRBs in the BATSE catalogs increases, we could expect factor of three improvements on 
the allowable population of Galactic GRBs. This should permit the number of ISCs to be 
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constrained. 


6. NEUTRON STAR VELOCITY DISTRIBUTION 


In the calculations above, we scaled the encounter rates to a fiducial NS space ve- 
locity u na = (100 km s -1 )vjoo, based on orbital calculations (Blaes & Madau 1993) of 
pulsars injected into the Galactic potential with a bimodal distribution of initial velocities 
ranging from 50 to 400 km -1 (Lyne, Anderson, & Salter 1982; Narayan & Ostriker 1990). 
However, new results indicate that pulsar initial velocities may be much higher. A recent 
re-evaluation (Lyne & Lorimer 1994) of the radio pulsar velocity distribution, taking into 
account a new distance scale, new proper motions, and selection effects, concludes that 
the mean pulsar space velocity at birth may be as high as 450 ± 90 km s -1 . There have 
also been suggestions (Kulkarni et al. 1994; Rothschild, Ivulkarni, & Lingenfelter 1994) 
that two SGR burst sources are associated with pulsars and neutron stars with extremely 
high space velocities. v n3 > 500 km s -1 . Therefore, it is worth considering the differences 
in our scenario for high-velocity NSs. 

Looking back to § 2 (eqs. [1] - [3]), we note that the NS - ESOC encounter rates 
will increase linearly with uioo, but the encounter durations will drop as Therefore, 

the expected total number of NSs intruding in ESOCs will remain invariant at around 
4000 f oc - However, the rate of comet - NS encounters in each of these sources (see eq. 
[7]) will decrease as If the SGR population turns out to be associated with NSs of 

abnormally large space velocities (u nj > 500 km s -1 ) their accretion rates from the ISM 
and their surface temperatures will drop. For accretion with v n3 c 3 , the mass accretion 

rate (see eq. [19]) and accretion luminosity scale as uj^g, and the surface temperature 
therefore will drop as v 10Q . The expected decrease in NS luminosity would result in 
considerably less radiative evaporation of comets during the shorter encounters (see eq. 
[13]). 

Thus, we conclude that, although the total number of “intruder NS” sources will 
remain the same throughout the Galaxy, the individual encounters between very fast, NSs 
and ESOCs will several times be less frequent, and the radiative ablation and subsequent 
accretion of cometary fragments will be much lower. It is therefore difficult to see how 
these encounters would be likely to make 7-ray bursts. As noted above, the more likely 
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sources, for our model of GRBs, are the older, slower population of NSs that have felt the 
gravitational potential of the Galaxv. 


7. SUMMARY 

The fact that the observed distribution of BATSE sources is both isotropic on the 
sky and volume limited (Meegan et al. 1992) has been used to suggest that most GRBs 
are not derived from Galactic sources (Paczynski 1991. 1993). Unless recent hints of small 
anisotropies in the distribution of intermediate count-rate GRBs (Quashnock & Lamb 
1993; Rutledge & Lewin 1993) turn out to be statistically significant (Schaefer 1993b), 
the apparent GRB isotropy argues that the majority of the bursters do not result from 
collisions between Galactic comets and NSs. 

We have shown that the rate of NS - comet collisions should be dominated by field NSs 
colliding with comets orbiting or ejected from other stars. Further, if collisions between 
comets and NSs produce GRBs, then most NS - comet collisions in the Galaxy should arise 
from collisions between field NSs and interstellar comets. These Extra-Solar Oort Cloud 
(ESOC) burst sources should produce space-correlated repeaters whose events exhibit a 
wide range of peak luminosities and luminosity time histories, suggestive of the soft gamma- 
ray repeaters (SGRs). Since the Galactic ESOC population is expected to reflect the 
geometry of the intermediate-population Galactic “thick disk” (Gilmore, Wvse, & Ivuijken 
1989), one expects a Galactic distribution for SGRs if comet - NS collisions are the burst 
mechanism. These general characteristics, and our rough estimate of 0.65 events yr -1 per 
source predicted for the NS - OC scenario, are qualitatively similar to the those of SGRs 
(Kouveliotou et al. 1993). 

Based on the expected event rates, we are able to constrain the fraction of stars with 
Oort Clouds like our own to be no more than 4%, unless either NS-comet collisions do not 
predominately produce 7 -rays or the 7 -rays are are highly beamed. If none of the detected 
SGRs known today lies within the BATSE-detectable distance scale for comet-NS events, 
the constraint on the fraction of stars with Oort Clouds like our own drops below 1 %. 
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